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First principles calculations show that H overlayers in- 
duce a honeycomb vacancy reconstruction on Al(lll). Above 
one ML H coverage, Al(lll) facets into surfaces vicinal to 
(111), with {100} faceted single and double steps and vacan- 
cies on the {111} terraces. These H induced reconstructions 
are stable because H binds more strongly to {100} than to 
{111} facets. Faceting of Al(lll) explains step bunching in 
H-mediated epitaxy of Al(lll). H-covered Al(100) is stable 
against faceting. 



Al is an important material for semiconductor devices. 
Alloyed with a few percent of Cu and Si, it is the domi- 
nant material for interconnects Q . Shrinking dimensions 
in chip design means that control of the Al deposition is 
more important. Al is usually deposited under conditions 
such that H is co-adsorbed on the growing film. However, 
H reduces the Al film's quality. Al(lll) grows layer-by- 
layer at 150 K in the absence of H ||. Calculations of 
Al surface diffusion barriers on flat and stepped Al(lll) 
predict smooth growth at even lower temperatures ||. 
With co-adsorbed H as a by-product of chemical vapor 
deposition (CVD), Hinch, Doak and Dubois report step 
bunching on Al(lll) Alms grown at 400-500 K @]. 

Since H is present in most commercial thin film deposi- 
tion processes, understanding the effects of H in epitaxy 
is of general interest. H does not only affect Al epitaxy. 
For example, H is an anti-surfactant in Si Upland Ni 
@ epitaxy, while on GaAs(100) J§ and C(100) |§ H im- 
proves epitaxial growth. 

Another major issue in H-solid interaction is the very 
common problem of H embrittlement. Materials that are 
usually ductile become brittle once they absorb H. Dur- 
ing H promoted fracture, H covers the newly formed sur- 
face. The energy of this H-covered surface can determine 
the fracture plane |l0| . 

To study the effect of H on the growth of, e. g., Al films, 
two surface properties are important: the mobility of Al 
on H-covered Al surfaces and the orientation dependence 
of the surface energy. A growing Al surface is rough if 
the lateral Al mass transport is hindered. A growing Al 
film facets if the surface mobility is high and if the to- 
tal surface energy can be reduced by forming facets off 
the growth direction 0. To determine surface mobility 
or surface energies the structure of the H-covered surface 
must be known first, because H can reconstruct surfaces 



]il[ . Indeed, H adsorption reconstructs Al(lll). As with 
Be(0001) [|nj, H induces a honeycomb vacancy recon- 
struction on Al(lll). In this honeycomb phase 1 ML of H 
atoms adsorbs on tilted bridge sites around the vacancies 
(see Fig. ^b). In this very open structure, Al diffusion is 
reduced. The effects of H modified surface diffusion will 
be studied elsewhere [[[2]. In this Letter we concentrate 
on understanding the step bunching observed in CVD of 
Al(lll) 0| . Our explanation is energetically favorable 
for Al(lll) to facet at high H coverage. The facet ori- 
entation is {311} or {211}, i.e. vicinal to Al(lll) with 
{100}-faceted steps (see Fig.[j]e+h). H induces vacancies 
and double steps on those facets. H-covered Al growing 
in the (100) direction does not facet, because H binds 
more strongly to Al(100) than to other Al surfaces. 

First-principles calculations The study of H adlay- 
ers on Al surfaces reported here is based on the local 
density approximation Jl3[] . Electronic wavefunctions are 
expanded in plane waves up to a cutoff energy of 12Ry. 
All atoms are represented by pseudopotentials Um. We 
model the Al surfaces by repeated slabs of 10-20 A thick- 
ness in orthorhombic or monoclinic supercells, separated 
by at least 6 A of vacuum. H is adsorbed on one side 
only. The adsorbed H atoms plus the top 4-6 A of the 
Al slab are relaxed. H adsorption energies are quoted 
relative to 1 Ry for free H. Special fc-points at a density 
equivalent to 500-5000 fc-points in the full Brillouin zone 
of a one atom cell are used. Details of the computational 
technique are published pjl lL 151. 

In earlier publications |3|,[ll] 16| we have accurately cal- 
culated properties of simple metal surfaces and of H- 
surface interaction. One example is the prediction of a 
novel type of H induced reconstruction, the honeycomb 
vacancy reconstruction of Be(0001) [jll). This reconstruc- 
tion is now confirmed quantitatively in experiment [fl7|| . 
We estimate that adsorption energy differences are sig- 
nificant if they are larger than 0.05 eV. Temperature and 
zero point motion effects are not included, however. Fur- 
thermore, it is uncertain how accurate calculated surface 
energies are, as surface energies are difficult to measure. 
Given these remaining uncertainties, we consider the suc- 
cess of our model in explaining a number of experimental 
findings for H on Al surfaces as the most important evi- 
dence for the validity of our theoretical results. 

H on flat Al(lll), Al(100) and Al(110) Properties 
of H adlayers on Al surfaces, like adsorption energy and 
coordination, depend on surface orientation. On the flat 
Al(lll) surface H prefers threefold fee sites with an ad- 
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sorption energy (E^) of 1.12 eV at 1/12 ML coverage (see 
Fig|l]a). The hep and bridge sites are 0.08 eV and the top 
sites are 0.22 eV higher in energy. That the fee site is the 
prefered threefold site for H on Al(lll) is significant and 
is discussed elsewhere At higher coverages the ad- 
sorption energy is about 0.05 eV higher (E ad = 1.17eV 
at 1/2ML and £ ad = 1.16eV at 1 ML). Thus, there is a 
weak attraction between H atoms adsorbed on Al(lll). 

On Al(100), H adsorbs preferentially at bridge sites 
(see Fig[l]f+g). The adsorption energy of H on Al(100) 
is higher than on Al(lll), with a maximum of 1.33 eV at 
1 ML. Below 1 ML, H adatoms attract each other (-E a d = 
1.22eV at 1/9 ML). Above 1 ML the H-H interaction is 
repulsive (£ ad = 1.26 eV at 1.5 ML and E ad = 1.09 eV 
at 2 ML). On Al(110), the H adsorption site is coverage 
dependent. Up to 1ML, H prefers the top site (E ad = 
1.19eVat 1/6ML, and E ad = 1.17eVat 1 ML). At 2 ML 
H covers the bridge sites of first and second layer (E ad — 
1.17eV). Thus the H-H interaction is weak. 

H at steps on Al(lll) It seems that adsorption of H 
on Al(lll), Al(100), and Al(110) is quite different and 
unrelated. However, a very close relation can be estab- 
lished by considering the adsorption of H at close packed 
steps on Al(lll). There are two types of close packed 
steps on fcc(lll) surfaces. One is a {lll}-microfacet 
(called B-step), where the step atoms have the same near- 
est neighbor environment as surface atoms on Al(110). 
The other step is a {100}-microfacet (called A-step). 

At B-steps H binds as on Al(110) at 1ML coverage 
and below jl8). The adsorption site is the top site, the 
adsorption energy, 1.19 eV, is nearly identical, and the 
H-H interaction energy is below 0.01 eV. At A-steps H 
binds at bridge sites. The H-H interaction is attractive 
with a maximum adsorption energy of 1.35 eV for one H 
adsorbate per step atom. Adsorption site and energy for 
H at A-steps is like for H on Al(100). 

The fact that H prefers the A-step over the B-step by 
as much as 0.16 eV has important consequences. One is 
that H adsorption on Al(lll) changes the equilibrium 
island shape. Without H the formation energy of A- 
steps is higher than that of B-stcps (0.248 eV compared 
to 0.232 eV per step atom) [||. Thus Al islands on clean 
Al(lll) are hexagonal, with longer B- than A-type edges. 
With enough H adsorbed to cover all step atoms, A-steps 
have less than half the formation energy of B-steps. This 
makes the equilibrium island shape triangular with only 
A-type edges. A change of island shape is not only of 
aesthetic interest, it also has been connected to a change 
in the growth morphology fll9| . The second consequence 
of the strong binding of H to A-steps is that H-covered 
Al(lll) transforms into a stepped surface with A-steps, 
as elaborated later. 

H adsorbed on Al adatoms Up to now we have dis- 
cussed the adsorption of H on flat and stepped Al sur- 
faces. Earlier studies of H adsorption on Be(0001) and 
Be(1010) show that H can cause vacancy reconstructions 



of a simple metal surface |11|,|12|,^6| . One experimental 
signature of the H induced reconstructions of Be(0001) 
is the high frequency H vibrational spectrum, which in- 
dicates strong bonds of the H to the reconstructed sub- 
strate jn|,|20|l. Similar modes are found for H on Al sur- 



faces 21 1 . However, the existing experimental results 



give no further hint concerning the structure of those 
possible reconstructions. To find a stable H-induced re- 
construction with theoretical tools, a large number of ge- 
ometries have to be studied. To speed up the search we 
identify stable H-Al units of which we then compose ge- 
ometries to study. This also gives insight into H-surface 
bonding transferable to other H surface systems. 

The driving force for H-induced reconstructions on 
simple metal surfaces is the desire of the H adsorbate 
to form strong (covalent) bonds to the substrate. These 
bonds are stronger if the metal surface atoms have a low 
coordination. Thus H adsorption lowers the energy of 
steps, adatoms and vacancies [0,^61. In fact, adatoms 
or vacancies might even become elements of the most 
stable surface geometry [ pd| . 

As on Be(0001), H binds strongly on top of Al adatoms 
on all three Al low index surfaces [|l6| . The adsorption en- 
ergy of H on top of a nearly isolated Al adatom is 1.64 eV 
on Al(lll), 1.74eV on Al(100) and 1.41 eV on Al(llO). 
Thus H is more stable on the adatom than in the most 
stable phase on the flat surface. The energy difference is 
0.47 eV on Al(lll), 0.41eV on Al(100), and 0.22eV on 
Al(110). On Al(lll) even a second H, adsorbed close to 
the adatom and on top of a neighboring surface atom, 
is stable. The total energy gain for 2 H atoms adsorbed 
at the adatom, compared to flat Al(lll), is 0.78 eV. The 
prefered binding of H to adatoms reduces the adatom 
formation energy, i.e. the energy to remove an Al atom 
from a bulk site and put it on the surface. Without H, 
this energy is 1.05eV on Al(lll), 0.56eV on Al(100), and 
0.29 eV on Al(110) at low coverage. With H, the forma- 
tion energies are reduced to about 0.27eV, 0.15eV, and 
0.07 eV. These values depend slightly on the H coverage. 
On Al(lll) the formation energy increases by 0.1 eV if 
the honeycomb phase (see Fig. [yb) is the H reservoir. 

Even though Al adatoms are more stable with H co- 
adsorbed than without H, we do not find any surface 
structure involving adatoms that has a lower energy than 
the flat surface at the same H coverage. However, the 
equilibrium Al+H concentration should be high on a 
H-covered Al surface. Assuming a room temperature 
Boltzmann distribution (kT = 0.025 eV), we estimate 
Al+H concentrations on H-covered Al(100) and Al(110) 
of 1/400 ML and above The Al+H ad-dimcr could 
be the precursor to Al hydride desorption. This would 
explain the observation that up to 10 % of the H adsorbed 
on Al(100) and Al(110) desorbs as Al hydride @. 

H-induced honeycomb reconstruction on Al(lll ) The 
most stable phase of H-covered Al(lll) is the \/3 x 
V3R30 honeycomb reconstruction. In the honeycomb 
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phase, every third Al(lll) surface atom is removed and 
H atoms cover the bridge sites connecting the remaining 
atoms (see Fig.|]b). The H adsorbates tilt by 26° in the 
direction of the closest fee site of the ideal surface. The 
H adsorption energy in the honeycomb phase {E^ d ) is 
1.21 eV per H, 0.05 eV higher than for 1ML of H on the 
flat surface. is defined as 

= -l/3(E^ - 3£ H - - 2El£) , (1) 

where E^ is the total energy per unit cell of a 5 + 2/3 
layer Al(lll) slab that has a honeycomb array of vacan- 
cies on one side. E n is the energy of a free H atom (1 Ry), 
i?5iay is the energy per unit cell of a 5 layer Al(lll) slab 
and -Egiay the energy of a 6 layer slab. This definition 
minimizes the errors introduced by finite slab thickness. 
H adsorption energies in other structures that involve re- 
moving Al atoms from the surface are defined similarly. 

What makes the honeycomb phase so stable ? The 
main reason is that all H atoms are coordinated like H 
atoms adsorbed at A-steps on Al(lll) (see Fig.|j]b+h). 
The H adsorption energy at A-steps is 1.35 eV, 0.19 eV 
higher than on flat Al(lll). In a simplistic model, the 
energy gain per V3 unit cell of having more stable 
H adsorbates would be 3x0.19eV= 0.57 eV. This more 
than compensates for the vacancy formation energy in a 
honeycomb array on Al(lll), which is 0.47 eV. 

Besides the honeycomb phase there are H-induced 
missing and added row reconstructions jll| at 1 ML H- 
coverage which are only 0.04 eV to 0.07 eV per surface 
atom higher in energy. These small energy differences 
and the small adatom formation energy should cause the 
H-covered Al(lll) surface to be disorded at room tem- 
perature. This could be the reason why only lxl pat- 
terns are detected in diffraction experiments. An alter- 
nate reason for the observed lxl patterns could be that 
the formation of the vacancy reconstructions is kineti- 
cally hindered. However, on Be(0001) the honeycomb 
structure already forms at 100 K Since Be(0001) 

and Al(lll) have similar surface diffusion barriers pjL6[|, 
we expect that the vacancy reconstructions of H-covered 
Al(lll) are activated at room temperature. 

The fact that H reconstructs Al(lll) and that it 
can even form a surface hydride at 2 ML coverage (see 
Fig. |l]c) might be the cause for the unusual H desorption 
from Al(lll): About half of the H desorbs as Al hydride 
and the H2 desorption rate is independent of coverage 
p3| . Details are presented elsewhere jl2| . 

H on surfaces vicinal to Al(lll) To find stable re- 
constructions of H-covered Al(lll) it is important to con- 
sider surfaces vicinal to Al(lll), especially those which 
have the energetically favorable A-steps. Furthermore, 
since H prefers (100) to (111), let us concentrate on vic- 
inal surfaces with small (111) terraces, i.e., Al(311) and 
Al(211). (Ill) terraces on Al(311) are two rows wide, on 
Al(211) they are three rows wide (see Fig. |l]c+li). 



As on Al(lll), it is favorable at high H coverages 
to form vacancies on the (111) terraces. These H sur- 
rounded vacancies are most stable if they are in the cen- 
ter of (111) terraces that are three rows wide. On ideal 
Al(311) the (111) terraces are too short, thus double steps 
are formed. The reconstructed (211) and (311) surfaces 
are very open and provide for a high concentration of 
favorable bridge sites for H adsorption. Nearly all H 
atoms are adsorbed on {100} facets or at A-steps. The 
H adsorption energy is between the values for the hon- 
eycomb reconstruction and the flat Al(100). The H cov- 
erage per area in those (311) and (211) phases, depicted 
in Fig.[l]e-|-h), is about 50% higher than ML coverage 
on Al(lll) or Al(100). Thus the (311) and (211) phases 
will be especially stable if the surface H concentration is 
high. Incidentially, the H saturation coverage measured 
on Al(lll), 0.19 A -2 pij, coincides with the coverage in 
the low coverage phase of the two (311) phases in Fig.[l]. 
We therefore suggest that in the measurement of the sat- 
uration H coverage on Al(lll) the surface was locally 
reconstructed with {311} and possibly {211} facets. 

Faceting of H-covered Al(lll ) The main goal of this 
Letter is to explain why H covered Al(lll) does not grow 
smoothly at 400 K but forms step bunches and is rough 
instead 0). We propose that these step bunches are 
{311}- and {211}-like facets which minimize the total 
surface energy of H covered Al(lll). To compare surface 
energies of H-covered surfaces we introduce a construct, 
the H chemical potential /xh • /iH is the energy of H in a 
fictitious infinite reservoir. Experimentally an increase in 
/in corresponds to an increase in H exposure or coverage. 
At large negative values of //h, H is expensive and no H 
adsorbs on the surface. Without H, the surface energy 
is £g Urf . If fJ-n is larger than the energy (— -E Had ) of H 
in the most stable phase (phase 1), then the surface is 
covered with a H coverage c\ and the surface energy is 

E snii (Mh) = ^surf ~ C l (MH + #H ad ) ■ ( 2 ) 

With increasing /in higher coverage phases (2, 3, . . . ) 
with smaller -E H ' ad " are formed. In principle transitions 
between different H phases can be used experimentally 
to correlate, e. g., H exposure and /xjj- 

Fig. ||a) shows the surface energies of five Al surfaces 
as a function of /^h- Without H, i.e. at small /uh, the 
Al(lll) surface is the most stable. With increasing /in, 
the Al(100) surface covered by 1 ML H becomes the most 
stable surface. At even higher /in, the reconstructed H- 
covered Al(311) and Al(211) surfaces take over. 

The surface energies plotted in Fig.|^a) are relevant to 
predicting the equilibrium crystal shape of H-covered Al. 
More importantly, these surface energies also determine 
the prefered crack plane in H induced brittle fracture, 
since the surface created in H induced brittle fracture is 
H-covered. We speculate that, as in Ga embrittlement 
of Al |To|], the prefered fracture plane is the (100) plane. 
However, further studies are necessary [12] . 
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Fig. ||b) shows the surface energies relevant for faceting 
of the Al(lll) surface. Except for Al(lll), all surface en- 
ergies are corrected for the extra surface area that would 
be created if the Al(lll) surface were to facet. For ex- 
ample, faceting into Al(100) creates 73% extra surface, 
which usually costs energy. Creating {311} or {211} 
facets increases the surface area only little (15% and 
6 %). Therefore, already at a H potential that slightly ex- 
ceeds that necessary to adsorb H on the Al(lll) surface 
(in the honeycomb phase), H-covered Al(lll) is unstable 
against forming trigonal pyramids with {311} or {211} 
faces. These pyramids might be identical with what HDD 
describe as step bunching at a growth temperature of 
400 K ||. Step bunching is less pronounced at 500 K j|. 
The major reason is that at 500 K the H desorption rate 
is much higher than at 400 K [fUS Thus the H con- 
centration on the surface and therefore /in, is reduced. 
Thus {311} and {211} faces are energetically less favor- 
able during CVD at 500 K than at 400 K. 

Conclusions In this Letter we give a consistent ex- 
planation why, in the presence of H, Al(lll) does not 
grow smoothly but is rough and forms facets Q. It is 
rough because H induces a honeycomb vacancy recon- 
struction with other reconstructions close in energy. Ad- 
ditionally, at higher H coverage, and if kinetics allow, 
H-covered Al(lll) transforms into reconstructed {311} 
or {211} facets because this minimizes the total surface 
energy. Al(100) is stable against faceting. This makes 
the (100) orientation a better candidate for smooth epi- 
taxial growth. 

The properties of the different H phases proposed here 
are compatible with or even suggest an explanation of 
measurements of the H saturation coverage, H2 and Al 
hydride desorption and the H vibrational frequencies. 
Clearly, the details have still to be filled in. Also, there 
are experimental hints that not only Al(lll) but also 
H-covered Al(100) and Al(110) reconstruct [glj, which 
deserves further investigation. 
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FIG. 1. H adsorption energy (eV) for H on five different 
Al surfaces as a function of coverage (1/A 2 ) and selected di- 
agrams of the atomic structure. 



FIG. 2. Surface energy (meV/A 2 ) of five Al surfaces as a 
function of the H chemical potential /j,h (eV). To the right 
the surface energies are corrected for the additional surface 
created if Al(lll) facets into a different surface. 
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